The Mesostructure of Polymer Collapse and Fractal Smoothing 
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We investigate the internal structure of a polymer during collapse from an expanded coil to a 
compact globule. Collapse is more probable in local regions of high curvature, so a smoothing of 
the fractal polymer structure occurs that proceeds systematically from the shortest to the longest 
length scales. A proposed universal scaling relationship is tested by comparison with Monte Carlo 
simulations. We speculate that the universal form applies to various fractal systems with local 
processes that promote smoothness over time. The results complement earlier work showing that 
on the macroscale polymer collapse proceeds by driven diffusion of the polymer ends. 
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Understanding the collapse of homopolymers from a 
flexible coil to a compact globule is a first step towards 
modeling the kinetics of molecular self-organization. It 
may be relevant to a description of DNA aggregation and 
the initial collapse of proteins from an expanded state to 
a molten globule from which the final ordered structure 
is formed. [1,2] We have performed scaling analysis and 
simulation of this transition to investigate kinetic effects 
during collapse. [3,4] The results suggest that the motion 
of the polymer ends plays an important role in kinet- 
ics because their motion is constrained only by a sin- 
gle bond. Along the contour monomers have two bonds, 
their motion is more constrained, and aggregation is more 
difficult. Thus, collapse for a long polymer occurs al- 
most as a one dimensional process where the polymer 
ends accumulate mass by moving along the contour of 
the polymer while accreting monomers and small aggre- 
gates. Encounters between monomers far apart along the 
contour to form rings are rare so they play no role in the 
collapse. As a result of the faster aggregation at the poly- 
mer ends the collapsing polymer on a macroscopic scale 
takes on a dumbbell like appearance. The few DNA fluo- 
rescence measurements that follow a single polymer col- 
lapse and its metastable states also indicate the special 
role of polymer ends. [5-7] However, a description of the 
internal structure of the polymer away from the polymer 
ends has not, thus far, been obtained. 

In this manuscript we consider the internal structure 
of the polymer during collapse, not including the ends. 
Our objective is to understand the local contour structure 
that consists of small aggregates and polymer segments 
between them. Our arguments generalize the consider- 
ation of the freedom of motion of monomers, because 
monomers found in straight segments are much more con- 
strained in their motion than monomers in curved seg- 
ments. This results in faster collapse in regions of high 



polymer curvature. 

We will focus on intermediate length scales between 
the size of the expanded polymer and the size of the col- 
lapsed aggregate. The length and time scales to which 
this analysis is relevant are between the size of the ini- 
tial coil, which scales as N v where N is the number of 
monomers and v = 0.6 is the Flory exponent, and the 
size of the final aggregate, which scales as N 1 / 3 (assum- 
ing a compact aggregate) . For long polymers these scales 
are well separated. During collapse, at these intermediate 
length scales, the internal structure of the final aggregate 
as well as of intermediate clusters that are formed should 
not be relevant. When convenient we can treat clusters 
as point objects, though this is not always necessary. The 
dynamic properties of cluster movement follows Stokes' 
Law-the diffusion constant of clusters decreases slowly 
with cluster size, D ~ R~ , where R is the radius of a 
cluster. This implies that the dynamics of clusters varies 
smoothly from that of the original monomers, and a uni- 
versal scaling behavior of the polymer during collapse 
should be found. By focusing on intermediate length 
scales, our results should be widely relevant to polymers 
with varied properties. While the eventual structure of 
the collapsed polymer depends in detail on monomer- 
monomer interactions, the separation of lengths scales 
implies that for a long enough polymer with a compact 
final aggregate, the details of these interactions should 
not be relevant to the kinetics of the collapse at early 
times. 

To characterize the collapse it is useful to compare the 
distance between two monomers with the contour length 
of the polymer connecting them. In conventional scal- 
ing the polymer end-to-end distance R is expressed as 
a function of the number of monomers N, or the num- 
ber of links in the chain L = N — 1. When aggregation 
occurs, the small aggregates that form, appearing like 
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beads on a chain, decrease the effective contour length of 
the polymer. We can define the effective contour length 
by counting the minimum number of monomer-monomer 
bonds that one must cross in order to travel the polymer 
from one end to the other. Bonds formed by aggrega- 
tion allow us to bypass parts of the usual polymer con- 
tour. Because we are not interested in the structure of 
aggregates we can neglect the difference between differ- 
ent kinds of bonds. In this way the effective number of 
links in the chain decreases over time. Thus, in order 
to study the internal polymer structure during collapse 
we investigate, via scaling arguments and simulations, 
the scaling of the end-to-end distance r(l,t) of internal 
polymer segments as a function of their effective contour 
length I. 

The equilibrium structure of the polymer before 
collapse — in good solvent conditions — is a self-avoiding 
random walk, where r ~ l v , and v = 0.6 in three dimen- 
sions. The contour length is proportional to the num- 
ber of monomers. During collapse, monomers are con- 
strained from aggregating with other monomers by their 
already existing bonds. A completely straight segment of 
polymer does not allow aggregation because no monomer 
can move to bond with another monomer. In contrast, 
highly curved regions are more flexible and monomers in 
these regions may aggregate. Aggregation in a curved 
region reduces the contour length and the polymer be- 
comes straighter, smoothing the rough fractal polymer 
structure. We therefore expect that the scaling exponent 
will increase over time. At long enough times the scaling 
will approach that of a straight line (r ~ I). However, 
this smoothing occurs first at the shortest length scales. 
In effect the polymer structure becomes consistent with 
a progressively longer persistence length. Assuming scal- 
ing, we anticipate that the polymer end-to-end distance 
for a polymer segment away from the ends of contour 
length I will follow the dynamic scaling formula: 

r = lf(t/n (1) 

The universal function f(x) is a constant for large values 
of its argument so that r ~ I (long times), and scales as 
x ( 1 - v )/ z f or small values of its argument, so that r ~ l v 
(short times). The short time regime described by Eq. (1) 
starts after an initial transient (a very short time regime) 
in which no new bonding has taken place. During the 
very short time regime the time dependence of the uni- 
versal scaling function does not apply. The usual scaling 
of the contour length and end-to-end distance persists 
until just after the very short time regime because the 
bonds that are formed initially do not form large rings 
and thus do not affect the large scale polymer structure. 
The short time regime begins with the first formation of 
individual bonds and lasts until the characteristic relax- 
ation time of the contour of length I. This time — the re- 
laxation time of the contour of length I — is the crossover 
time between the short and long time regimes which fol- 
lows a scaling law r ~ l z . The dynamic exponent z is 



assumed to be consistent with conventional Zimm relax- 
ation, z — 2>v . Finally, we can also rewrite this scaling 
relation in terms of the number of monomers n in a poly- 
mer segment. Since the average mass along the contour 
is M ~ n/l and M follows power law scaling[3] M ~ t s — 
we substitute I ~ nt~ s in Eq. (1) to obtain r(n, i). 

We emphasize that kinetic effects become important 
for collapse of polymers in poor solvent, after equilibra- 
tion in good solventthe result of a quench in solvent affin- 
ity or temperature below the thermodynamic transition 
at 9-solvent conditions. Close to the 6-point a mean 
field argument where kinetics do not play a significant 
role is likely to apply. [8,3] In contrast, we will approx- 
imate the collapse by a completely irreversible model 
where no disaggregation occurs. Because the scaling vari- 
able that determines the effective distance from the 0- 
point is A 1 / 2 AT, long lengths are equivalent to small 
temperatures, and microscopic reversibility becomes ir- 
relevant at long enough length scales. We are thus con- 
sistently adopting a description that is valid for lengths 
longer than the microscopic regime. We further restrict 
our study to diffusive monomer motion and short-range 
interactions. 

The scaling relationship, Eq. (1), was tested by Monte 
Carlo simulations. These simulations in part include 
the effects of hydrodynamics during collapse by scaling 
the diffusion constant of aggregates according to Stokes' 
Law. The simulations are based on the two-space lat- 
tice Monte Carlo algorithm [3, 4, 9, 10] developed for sim- 
ulating high-molecular-weight polymers, and shown to 
be significantly faster than previous state-of-the-art tech- 
niques. [9, 10] In the two-space algorithm odd monomers 
and even monomers of a polymer are distinct and may 
most easily be described as residing in two separate 
spaces. Each monomer occupies one cell of a square 
lattice. Both connectivity of the polymer and excluded 
volume are imposed by requiring that, in the opposite 
space, only the nearest neighbors along the contour re- 
side in the 3x3x3 neighborhood of cells around each 
monomer. Motion of monomers is performed by Monte- 
Carlo steps that satisfy the polymer constraints. Since 
adjacent monomers (and only adjacent monomers) may 
lie on-top of each other, the local motion of the polymer 
is flexible. Despite the unusual local polymer properties 
the behavior of long polymers is found to agree with con- 
ventional scaling results. 

The polymer is initially relaxed into an equilibrium 
configuration using a fast non-local "reptation" Monte 
Carlo algorithm. Monomers are randomly moved from 
one end of the polymer to the other, which, for equilib- 
rium geometries, provides equivalent results to the local 
two-space dynamics. 

Collapse of the polymer is then simulated using lo- 
cal diffusive Monte Carlo dynamics, but without the ex- 
cluded volume constraint. Simulations of a variety of 
models indicate that excluded volume does not signif- 
icantly affect the kinetics of collapse. [4] Monomers are 
no longer stopped from entering the neighborhoods of 
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FIG. 1. 'Snapshots' of the collapse of a single 

homo-polymer of length N = 500 monomers in 2-d using the 
two-space algorithm. The plot is constructed by placing dots 
of area M 1//2 for an aggregate of mass M. This does not reflect 
the excluded volume of the aggregates, which is zero during 
this collapse simulation. The frames demonstrate the pro- 
cess of local smoothing that occurs progressively from short 
to long length scales. The behavior of the ends is discussed 
in Refs. 3 and 4. The primary effects of hydrodynamics are 
included in the simulations by applying Stokes' law to the 
diffusion of aggregates. 



other monomers; they continue to be required not to 
leave any neighbors behind. This allows monomers in 
the same space (odd or even) to move on top of each 
other, and thus aggregate. Aggregates of any mass oc- 
cupy only a single lattice site, and are moved as a unit 
by the same dynamics used for monomers. The mass 
of an aggregate is set equal to the number of monomers 
located at that site. The probability of moving an ag- 
gregate is adjusted to be consistent with a diffusion con- 
stant that scales by Stokes' law for spherical bodies in 
3-dimensions, D ~ M -1 / 3 . This represents the effects 
of hydrodynamics on individual aggregates, but does not 
include coupling of motion of different aggregates. One 
time interval consists of attempting a number of aggre- 
gate moves equal to the number of remaining aggregates. 

The end-to-end distances of polymer segments, r, were 
measured as a function of the effective contour length, I. 
The effective contour length is the minimum number of 
links along the polymer that connect a monomer at one 
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FIG. 2. Plot of the internal polymer segment end-to-end 
distance r, as a function of segment contour length I, and 
collapse time, t, in both (a) two dimensions (2-d), for t =0, 
25, 50, 100, 200, 400, 800, 1600 and (b) three dimensions 
(3-d), with t =0, 25, 50, 100, 200, 400, 800. In both cases the 
polymer contained 500 monomers and results were averaged 
over 200 collapses. 



end of the segment with the other end of the segment. 
Since an aggregate occupies only a single lattice site, in- 
terior bonds of the aggregate need not be counted, and 
it can be treated like a single monomer. The end-to-end 
distance is not exactly the Euclidean distance between 
the ends. Instead it is correctly defined as the mini- 
mum number of links needed to connect the two ends 
by any curve in space. Due to the underlying lattice in 
our algorithm a Manhattan metric, with the inclusion of 
diagonals, is appropriate. 

In addition to simulations of polymer collapse in three 
dimensions we also performed simulations of polymers 
whose motion is confined to two dimensions, which are 
convenient for pictorial illustration (Fig. 1). These are 
not conventional two dimensional simulations because of 
the well known problems with hydrodynamics in two di- 
mensions. [11] Instead, they represent the dynamics of a 
polymer confined at an interface (e.g. between two flu- 
ids). Thus the polymer is confined to two dimensions 
while the hydrodynamics is three dimensional. In this 
case we have v = 0.75, z = Zv and D ~ M -1 / 2 . The 
frames in Fig. 1 illustrate contour smoothing. Starting 
with short length scales, the polymer becomes progres- 
sively smoother and approaches a straight line. 

In Fig. 2 we show log-log plots of contour length (I) 
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FIG. 3. Plot of the scaling exponent d In r/d In I as a func- 
tion of time for different segment contour lengths /. This 
figure was obtained from Fig. 2 by calculating average slopes 
over a segment length of 20. 

versus end-to-end distance (r) for both 2-d and 3-d simu- 
lations. Initially the results are consistent with r ~ l v for 
a self avoiding random walk. As time progresses the poly- 
mer becomes smooth resulting in a slope that approaches 
1. The asymptotic behavior can be seen to occur earlier 
at the shortest length scales. 

Fig. 3 shows the derivative, obtained from finite differ- 
ences, as a function of time for different segment contour 
lengths. For all segment lengths the derivative starts at 
approximately v and approaches 1 as the collapse pro- 
ceeds. The rate of collapse becomes progressively slower 
as I increases. The scaling relation, Eq. (1), predicts that 
the relaxation time will scale with I as l z . Fig. 4 shows 
the data following rescaling. r/l is plotted against the 
rescaled time, t/l z . The generally good coincidence of 
the different curves confirms that the simulation obeys 
Eq. (1).[12] An attempt to use an asymptotic scaling ex- 
ponent, r <~ l u , with u = 0.95, led to a visibly poorer fit, 
as did small variations in the exponent z. 

The excellent agreement with the expected universal 
scaling relationship using the hydrodynamic exponent, 
z = "Sv , may be fortuitous because the simulations do 
not contain the full effects of hydrodynamics. Specifically 
they contain only the effect of hydrodynamics on individ- 
ual aggregates and not the coupling between aggregate 
motion. 

In summary, we have found that polymer collapse dis- 
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FIG. 4. Plot of the rescaled end-to-end polymer segment 
distance, r/l, as a function of the rescaled time, t/l z . Data 
for both 2-d and 3-d are shown. The coincidence of the curves 
is consistent with validity of the universal scaling relationship, 
Eq. (1). 

plays a process of fractal smoothing that occurs first at 
the shortest length scales. Our simulations were found 
to be in good agreement with a universal scaling rela- 
tionship. It is interesting to speculate that this may also 
apply to other fractal systems where local smoothing pro- 
cesses occur. Several groups have attempted to measure 
the self-affine scaling of horizontal transects of mountain 
ranges. [13,14] They found that a unique fractal dimen- 
sion, Dh cannot be assigned, but that the effective frac- 
tal dimension decreases with length scale. For example 
Dietler and Zhang [13] have performed calculations for 
Switzerland, an area of 7 x 10 4 km 2 , with a resolution 
of 100 m. They obtained Dh ~ 1-43 at length scales 
below approximately 5 km, and Dh ~ 1-73 for larger 
length scales. The data points could also lie on a con- 
tinuous curve rather than two distinct scaling regimes. 
Thus the landscape appears smoother at shorter length 
scales. Short range smoothing may arise from processes, 
such as weathering, that also give rise to short range cor- 
relations. Various fractal biological systems formed as a 
result of an initial developmental process may also suffer 
smoothing as part of aging. 

Since this work[15] was completed, a number of other 
works have explored the kinetics of collapse using simula- 
tions, analytic treatments and scaling arguments. Tim- 
oshenko, Kuznetsov and Dawson[16] studied the kinetics 
of collapse using Monte Carlo simulations and a mean- 
field "Gaussian self-consistent" approach. Their Monte 
Carlo simulations are based upon an underlying lattice 
model which is similar to ours, however, they do not move 
aggregates as a unit. Since they only move individual 
monomers, Monte-Carlo rejection of moves causes the 
diffusion constant of aggregates to decreases very rapidly 
with aggregate size (naively, it decreases exponentially 
[17]). By contrast, in a fluid, collective motion results in 
Stokes' law diffusion, which we have included in our sim- 
ulations. The slow diffusion of aggregates in their simu- 
lations cause their results to be distinct from ours. From 
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their figures it appears clear that aggregates tend to pin 
the polymer contour. Their Gaussian self-consistent ap- 
proach is analytically elaborate, however, it is not clear 
from their analysis whether it treats correctly the diffu- 
sion of clusters. Moreover, since some equilibrium scaling 
laws are not correct in this method it is hard to evaluate 
whether the kinetic properties are correct and their anal- 
ysis does not clarify this point aside from the claim that 
the analytic results are in aggrement with their Monte 
Carlo simulations. 

Buguin, Brochard-Wyart and de Gennes[17] have pre- 
sented scaling arguments based on a model of local clus- 
ters "pearls" forming during collapse close to the 0-point 
in the mean field regime where surface tension is the driv- 
ing force of collapse. Pitard[18] has further considered 
the dynamics of collapse in this mean field regime by dis- 
cussing the effect of tension along a polymer contour be- 
tween two clusters (pearls) and extended the arguments 
to considerations of a string of clusters. These papers 
refer to a different regime (i.e. the mean field regime) 
than our analysis. Within this regime they provide com- 
plementary insights about the structure of clusters or 
pearls during collapse and the formation of a globule, 
which is important both to the kinetics of collapse and 
to the eventual structure of the aggregate that is formed 
at the end. It is worth noting that our simulations do 
not allow monomer motion along the polymer contour 
which can allow monomers to leave and join aggregates. 
The distribution of cluster sizes may be affected by such 
motion. We note, however, that the essential results of 
this paper should not be changed by redistribution of 
monomers along the contour, and resulting change in the 
distribution of the sizes of clusters, because they only 
affect the distribution of diffusion constants which vary 
only weakly with aggregate size. Moreover, the scaling 
law Eq. (1) does not refer to aggregate size and should 
not be affected. 

Finally, Kantor and Kardar[19] have investigated the 
properties of charged polymers and find their compact 
form exhibits a necklace shape with end aggregates and 
intermediate aggregates forming as a function of the 
charge density. These results also display some interest- 
ing similarities to collapse behavior and further research 
may reveal a connection between their results and the 
studies of collapse. 

We would like to thank A. Grosberg and M. Kardar 
for helpful discussions. A referee is to be acknowledged 
for pointing out that the universal scaling law does not 
apply in the limit t — > before the first bonding events. 
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